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FIG. 1: Sketch of the setup: M stands for mirrors and B for 
balanced beam-splitters. A non-degenerate parametric ampli- 
fier (down-converter DC) generates two-mode squeezed vac- 
uum in modes ài and òi . These are fed into an interferometer, 
with a phase shifter $; its output modes 03 and 63 are de- 
tected by homodyne detectors H a and H;, whose signals are 
multiplied. Note, that neither the amplifier's pump-laser nor 
the làser for homodyne detection is shown, also, typically, the 
polarization of modes ai and òi does not conform and has to 
be adjusted. 



The use of balanced homodyne detection removes the 
detection problems mentioned above. Because only well 
established technology is required |Ï8|, |Ï9|, a proof-of- 
principle experiment will be immediately possible. 

In order to derive our main result ([!]) we follow the con- 
ventions of reference |0J : In the Heisenberg picture the 
action of the parametric amplifier is described by pho- 
ton operator transformations à\ = Uào + Vb\ and b\ = 
Ub + Và\ where U — coshG and V = — i exp(i£) sinhG 
with the single pass gain G — g\E p \L and a relative phase 
£ which we will assume to be zero. L is the interaction 
path length, E p the pump laser's amplitude, and g the 
gain coeficient proportional to the nonlinear suscepti- 
bility °f the down-conversion medium DC. Beam 
splitter Bi is described by à 2 = exp(ï$)(ài — ib±)/y/2 and 
62 = {—ià\ +òi)/v / 2; note that the interferometric phase 
shift <ï> in arm hi is included. The action of the beam 
mixer B 2 is analogously described by a 3 = — xb·^j·spï 
and 63 = (— ià2 + &2)/v2 and the total transformation 
thus reads 



«3 



l -^-{Uà + Vbl) + l -^—(Ub + Val) (2) 



-{Uào + Vbl) + ~^^(Ubo + Và\) .(3) 



2i 

Since we assume that modes ào and bo are in the vacuum 
state, two mode squeezed vacuum in modes Si and b\ re- 
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FIG. 2: Signal (P) and the square root of the signals second 
moment (P 2 ), see eqns. (^|) and (^), give us an idea of the 
behavior of the noise when measuring the homodyne current 
product (fel). 



sults, parameterized by the squeezing or gain parameter 
G. The corresponding intensity (N) is M 



(N) = (0% + blb 3 ) = (ajài + b\h) = 2sinh(G) 2 . (4) 

It is well know that balanced homodyne detection mea- 
sures the quadrature components of the monitored fields. 
We assume a relative phase of zero between local oscilla- 
tor and our interferometric modes à 3 and 63 . In this case 
the photo currents of detectors A and B are proportional 
to the expectation vàlues of à 3 + 03 and Ò3 + 63 § ü- 
The product P of the photo-currents is the signal we are 
interested in, it amounts to 



(P) 



((a 3 + à 3 )(bl + 63)) = sinh(G) cosh(G) sin(2$).(5) 



Note, that we observe a double period in the phase in- 
terval $ = 0, 2w in Eq. (||) and in FIG. || because our 
signal stems from the product of two homodyne currents. 
The corresponding second moment (P 2 ) is 



(P 2 



l + [ 



\ + cos(2$) - \ cos(4$)] (^L + (JV)) ,(6) 

where we used the intensity expression d). This yields 
the Standard deviation a = J (P(í>) 2 ) - (P(<f>)) 2 

= ^(^ + W )[| + cos( 24 )-^Í)] + 1 ,(T) 

which is minimal for ^ m in = 4> = t/2. Consequently 
the associated Standard expression for the phase resolu- 
tion limit 5<j> = CT($)/|aP/9$|$ =ff/2 is 1/|<9P/<9$|<ï, =0 = 
í/y/(N) 2 + 2(N) ps l/(iV). This result seems to indi- 
cate that we can reach the Heisenberg limit since the 
minimal detected phase difference 5(j) ps 1/ (-/V) . But an 
inspection of the behavior of the second moment of the 
signal in FIG. |2| shows that the noise varies greatly in 
the vicinity of the optimal point <fi = ir/2. We therefore 
have to analyze the behavior of the noise-valley around 
<fi more closely. It turns out that the ràpid growth of 
noise away from the optimal point does not let us achieve 



FIG. 3: Logarithmic plot of minimum phase spread log 10 (5w) 
as a function of the single pass gain G determined from (H) 
(curved black line) and of four times the inverse intensity 
log 10 (4/(2 sinh(G) 2 )) (straight red line) confirming that óip = 
4/(2sinh(G) 2 ) =i/{N). 



Heisenberg-limit resolution but the gradient of the slopes 
is sumciently low to allow for a reduced phase resolution 
that s cales like the Heisenberg limit, namely, according to 
our main result d) . Note, that a similar problem was en- 
countered in referenec |3j which was resolved by the stip- 
ulation that the interferometer acted 'phase-conjugated', 
meaning, when arm ai lengthens 62 contracts by the 
same amount. In the present case this solution does 
not help and we have to accept a diminished perfor- 
mance. To derive our limit ([!]), let us remind ourselves 
of the Standard derivation for the noise-induced phase- 
spread that limits interferometric resolution. Assuming 
that we encounter a noisy signal P( < í > ) ± with Stan- 
dard deviation a we want to be able to teli the param- 
eter <f> apart from cf> + Ac/>. We therefore require (as- 
suming, for definiteness, that P > and growing with 
increasing <E>) that, according to the Raylcigh-criterion, 
P(4>) + ^f- < P{(f> + A(j)) - g (»+ A ») , Approximating 
P(<p + Acj)) « P(<£) + A<t> ■ dP((j))/d(j), assuming equal- 
ity of left and right hand side in order to determino 
the smallest permissible 54> and that the variance does 
not change appreciably a((f> + Acj)) sa <j((f)) this yields 



the Standard expression for the phase resolution limit 
S(f) = (j(4>)/\dP/ 8^1^=4,. In our case, however, we need 
to look at an expression which accounts for the chang- 
ing variance; we therefore have to include both variances 
<j(ip) and a(tp + 5(p), according to the above discussion 
this leads to the modified criterion 



5cp = 



a(ip) + a((p + dtp) 



1 



\dP/d^-_ 



(8) 



Choosing the optimal working point ip — n/2, this yields 
an implicit equation for 5<p which is not too easy to solve 
in the general case but for sufncicntly high intensitics 
(G > 2.5 <-> (N) > 2sinh(2.5) 2 ps 73 photons) we find 
S(p = A/(N). This is illustrated by FIG. | and can be 
verified by direct substitution into (||). Because in our 
scheme the noise is phase sensitive it only works at par- 
ticular phase settings (odd múltiples of 7r/2, see FIG. ||) 
and our setup has to include a feedback mechanism - not 
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FIG. 4: Logarithmic plot of minimum phase spread \og 10 (Sip) 
determined from ([]) (curved black lines) and of four times 
the inverse intensity log 10 (4/(2 sinh(G) 2 )) (straight red lines) 
comparing various loss mechanisms: a) ai = /3i = 7r/300, b) 
«2 = f3 2 = tt/300, c) Ai = 1/300000, d) A 2 = -1/3. Note the 
different orders of magnitude in the vàlues of Ai and A2 and 
the fact that 'd)', surprisingly, shows a slight improvement in 
performance for a large negative imbalance A2 = —1/3. 



mentioned in FIG. |]. 

Robustness and further increase in sensitivity: 
Having shown that our scheme allows for Heisenberg- 
limit-like scaling in interferometric sensitivity we would 
also like to look at its sensitivity to experimental imper- 
fcctions. Balanced homodyne detection amplifies quan- 




For strong fields 
lÏ9l E2| and will 



tum features to the classical level 
detector losses can be kept small 
therefore not be discussed further. 

More importantly, losses and imbalanccs in the 
state preparation and interferometric part of the setup 
sketched in FIG. [ï] deserve consideration. The main ques- 
tion we want to address is whether the introduction of 
experimental imperfections leads to a gradual loss of per- 
formance or whether we might be unlucky and a quali- 
tative change in behavior results from any minute im- 
perfection. It turns out that the former is the case, yet, 
experimental demands on the state preparation part of 
the setup are very high. 

FIG. | compares the various cases for losses and im- 
balances and shows that the system is more forgiving for 
losses in the interferometer part than in the state prepa- 
ration part: the utilized quantum state has to be pre- 
pared with great skill but the scheme is comparatively 
robust to imperfections of the interferometer. When all 
imperfections are studied simultaneously their effects add 
up, i.e., tend to be dominated by the largest effect(s). 

Let us first consider losses in the state preparation 
part of the setup, i.e. losses in modes ài and bi extend- 
ing from inside the crystal to the first beam-splitter Bi. 
They are described by the mode-transformations à\ 1— > 
cos(ai)ài + sm(ai)òi and b\ 1— ► cos(/?i)fei + waifi\)V\ 
plus subsequent tracing over the loss modes (not men- 
tioned) and the admixed vacuum modes U\ and V\. It 
turns out that the qualitative picture does not depend 



much on the details such as whether the loss parameters 
ol\ and (3i are equal or the losses oceur in one channcl 
only. Thus, with experiments in mind, let us assume 
symmetric losses, namely a± = f3\ = 7r/300 w 0.01 lead- 
ing to 0.0001 = 0.01 % losses in both channels. For 
example, mode-mismatch at the beam splitter Bi leads 
to such symmetrical admixture of vacuum. This scheme 
is very sensitive to losses in the state preparation part 
of the setup and shows saturation of performance, see 
FIG. D a): to gain an order of magnitude in performance 
«i and (3\ have to be decreased by half an order of mag- 
nitude, namely, 5ip saturates at about 9af . 

Losses in the interferometric part of the setup (modes 
a-i and 62) are analogously described by loss parameters 
«2 and P2 which parameterize the admixture of two more 
vacuum modes u% and £2 to the path modes ài and 62- 
FIG. ^ b) illustrates the greater tolerance of our scheme 
to losses in the interferometer part of the setup. For 
the same loss vàlues, a.2 = /3 2 = tt/300 as for ol\ and (3\ 
before, the scheme shows a relatively better performance, 
indeed, the performance does not show saturation at all. 
Instead, at the threshold intensity » 4/ (9a 2 ) it switches 
from 5<p « 4/(iV) to the poorer scaling 5<p = k/ yj {N), 
maintaining the ground it has gained. Namely, beyond 
the threshold intensity we find 5ip w 6 \ot2\/ \J (N). 

Analogously to the case of losses, the scheme is also 
much more sensitive to imbalances of the first beam- 
splitter Bi than to those of B2, which is described by 



B 2 







«2 









-cos(f + A 2 ) isin(f + A 2 ) 
-<BÍn(|+A 2 ) cos(f + A 2 ) 



(9) 



conforming with the case A2 = used in the derivation 
of Eqs. (0) and (||). The imbalance in transformation 
Bi is described, in full analogy, by an imbalance angle 
Ai. Similarly to the case of a\ a non-zero Ai leads to 
saturation: for positive imbalances the saturation level is 
Sep i=y 4Ai and it is 6ip ~ 12|Ai| for negative Ai. 

It turns out that variation of A2 modifics the cocffi- 
cient k but not the scaling exponent in 5tp — n/(N); as 
mentioned above, for A2 = we find 5tp = 4/ (N). Note 
that, surprisingly, a large (negative) imbalance A2, as is 
displayed in FIG. |^ d), yields a small increase in perfor- 
mance quality (k is being reduced) . I cannot explain this 
finding and I think it deserves further investigation and 
might even lead to a trick to reduce the scaling reported 
here down to the Heisenberg-limit Síp = 1/(N). Varia- 
tion of the value of the imbalance parameter A2 alone, 
leads to an optimal value for the imbalance of approxi- 
mately A2 « —0.2375= — 13.61° and to our central result 
Eq. (|l]). This probably is the best our scheme can offer. 

Over recent years a consensus has emerged that a 
sharp photon number distribution is needed to reach 
the Heisenberg-limit || ||, || It was there- 

fore even concluded that the perceived need of a sub- 
poissonian photon number distribution renders the two- 
mode squeezed vacuum state unsuitable for interferome- 
try because of its super-poissonian thermal photon num- 
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ber distribution O, [ïl| ; in the light of these claims our 
central result Eq. (pj) is rather surprising. 

Note, that we did not discuss a criterion for the 
power of the pump beam driving the parametric down- 
conversion source and of the power needed for the strong 
local oscillator fields necessary to perform the balanced 
homodyne measurements H a and H;,, see FIG.§ Ifthis is 
included, the effective performance of our scheme could 
be reclassified as less efficient, yet, it remains a scheme 
with Heisenberg-limit-like scaling. The penalty to pay is 
not too large for large intensities because the local oscilla- 
tor's shot-noise-to-signal-ratio diminishes with increasing 
signal strength thus yielding very accurate homodyning 
signals. In this context I would also like to mention that 
there are promising recent ideas for efficient and bright 
down-conversion sources p3| . 

Also note, that the considerations of this paper might 
turn out to be of importance for atom-beam interferome- 
try |Ï5|| since four-wave mixing has been reported yielding 
correlated atom beams in states similar to the two-mode 
squeezed vacuum states discussed here Q . 



Conclusions: We have found that bosonic two-mode 
squeezed states can be used in an interferometcr to 
achieve phase resolution near the Heisenberg-limit, see 
Eq. (0). This only works at particular phase settings, the 
noise is phase sensitive and the setup therefore needs a 
feedback mechanism. The degrading influence of exper- 
imental imperfections is analyzed and it is shown that 
requirements on the state-preparation part of the setup 
are very stringent. On the other hand our scheme is more 
robust with respect to imperfections of the interferome- 
ter part of the setup and it does not suffer from single- 
photon detection problems because it relies on balanced 
homodyne-detection. 
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